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Wo present a study of 649±35 e^e~— >cc events produced at 10.6 GcV containing both a ylj 
baryon and a A'^ antibaryon. The number observed is roughly four times that expected if the leading 
charmed hadron types are uncorrelated, confirming an observation by the CLEO Collaboration. We 
find a 2-jet topology in these events but very few additional baryons, demonstrating that the primary 
c and c are predominantly contained in a correlated baryon-antibaryon system. In addition to the 
charmed baryons we observe on average 2.6±0.2 charged intermediate mesons, predominantly pions, 
carrying 65% of the remaining energy. 

PACS numbers; 13.66.Bc, 13.87.Fh, 13.60.Rj 



Baryon production in high-energy jets from e+e an- nihilations has presented a series of challenges to our un- 
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derstanding of strong interactions. Its observation led to 
the competing notions of 'primary' and 'local' baryon 
correlations In the former, the e+ and e" anni- 

hilate into a primary diquark-antidiquark, rather than 
a quark-antiquark, pair. The diquark and antidiquark 
then hadronize into jets containing a leading baryon A^i 
and a leading antibaryon 7V2, respectively, but no other 
(anti)baryons. Ni and A^2 would then share two quark 
flavors and typically have high, antiparallel momenta and 
large values of variables characterizing their separation, 
such as invariant mass or rapidity difference |Aj/|, where 
y = 0.5 ln[(£J + — E is the baryon energy, 

and p|| is the projection of its momentum on the thrust 
axis. Alternatively, an N1N2 pair might be produced lo- 
cally, in an individual step of a hadronization cascade, 
with a smaller value of |Ay|. Most experimental studies 
of baryon-antibaryon pairs have shown \Ay\ distributions 
that peak at small values 0. 

Several mechanisms to describe baryon production 
and correlations have been implemented in Monte Carlo 
hadronization models [1]. In the JETSET 4] color- 
flux-tube model, a tube break can result in a diquark- 
antidiquark (rather than qq) pair, producing an N1N2 
pair locally. An intermediate meson is introduced be- 
tween A^i and A^2 with some probability (50% by de- 
fault [sl) to match the measured \Ay\ distributions. In 
the HERWIG Q model, an individual, color-singlet clus- 
ter may fragment into a baryon-antibaryon pair but not 
a multi-body state with additional mesons. The model 
does not reproduce the measured \Ay\ distributions when 
tuned to other observables [Ij. The UCLA area-law 
model includes N1N2 pairs with any number of interme- 
diate mesons, and suppresses higher-mass intermediate 
meson systems by means of a tunable parameter. 

Direct evidence of primary production and/or inter- 
mediate mesons would be of great interest, but previous 
searches for the latter using three-particle correlations \^ 
or baryon flavor correlations were generally inconclu- 
sive. 

At center-of-mass (cm.) energies ^/s much larger than 
four baryon masses, the assumption of local baryon num- 
ber conservation implies that an e+e"— )■ qq event contain- 
ing a leading baryon A'^i in the 9 jet and a leading an- 
tibaryon N2 in the q jet must also contain an antibaryon 
A^3 in the q jet and a baryon N4 in the q jet. However, 
if the N1N3N4N2 mass is a large fraction of ^/s, these 
four-baryon events would be suppressed and other pro- 
cesses might be visible — in particular, primary baryon 
production events with exactly two baryons, one in each 
jet. At i/s « 10 GeV, charmed (c) baryons are of par- 
ticular interest, since any high-momentum c or c baryon 
must be a leading particle in an e"'"e~— J-cc event, and any 
NciN3NiNc2 mass exceeds 6.5 GeV/c^. The CLEO Col- 
laboration reported an excess by a factor of 3.5 ± 0.6 
in the number of events at ^/s = 10.6 GeV with both a 
A'^ and a A~ , where their expectation is derived assum- 



ing local baryon number conservation in the JETSET 
model and from observed events with a A'^ and a D~ 
or D meson. This excess is evidence that the baryon 
production is correlated between the c and c jets and 
is consistent with primary baryon production, but does 
not exclude the possibility of local baryon production 
with correlation between the jets. The two cases can be 
distinguished experimentally: local production would re- 
quire an additional baryon and antibaryon (A'4 and N3) 
in the event, so events with exactly one ylji", exactly one 
A~ , and no additional baryons would imply primary pro- 
duction. CLEO investigated this and did not observe a 
strong signal for additional protons in the A'^A'^ candi- 
date events, but due to a limited data sample and the 
lack of a limit on additional neutrons they were unable 
to exclude local baryon production. 

In this paper we exploit the particle identification ca- 
pabilities of the BABAR detector jTl[ to select a sample 
of A'^A~X events in which the yl+ and A~ are pro- 
duced at high momentum in opposite hemispheres, and 
study their characteristics in detail. We use 220 fb"^ of 
data collected at ^/s= 10.54-10.58 GeV. We identify the 
charged tracks in the X system, looking for additional 
(anti)protons, and search for higher-mass baryons that 
could be a source of the A^A~X events. We consider 
charged tracks measured in the silicon vertex tracker 
(SVT) and drift chamber (DCH), and identified as pi- 
ous, kaons or protons using the DCH and the detector of 
internally reflected Cherenkov light. The identification 
algorithm used here [13, 113] is over 99% efficient for pions 
and kaons (protons) within the acceptance with momenta 
between 0.15 and 0.5 (1.2) GeV/c, with misidentification 
rates below 0.5%. At higher momenta it remains over 
90% efficient, with misidentification rates generally be- 
low 1%. 

We construct yl+ candidates in the pK ~Tr~^ and pK'^ 
decay modes and A~ in the corresponding charge- 
conjugate modes. We consider a pair of oppositely 
charged tracks as a K'^ — >■ tt+tt" candidate if a vertex 
fit returns a with a confidence level (CL) exceeding 
0.01, the vertex is displaced by 2.5-60 cm from the inter- 
action point (IP) calculated for each event from the set of 
well-measured tracks in the SVT, the angle between 
the Kg candidate's momentum and the IP-to-vertex di- 
rection satisfies cosOks > 0.97, and the tt+tt" invariant 
mass is in the range 491.8-503.8 MeV/c^. All combina- 
tions of a ifg and a well-measured (>15 hits in the DCH 
and >5 in the SVT) proton are considered A'^ — )■ pK*^ 
candidates. A combination of well-measured p, K~ and 
7r+ tracks is considered a A'^ ^ pK^ir^ candidate if its 
vertex fit yields CL>0.001. 

We require p*, the momentum of the A'^ candidate 
in the e+e^ cm. frame, to exceed 2.3 GeV/c, so that 
the rate of A'^ from T(4S') decays [13, [l3| is negligible. 
We select events containing at least one yl+ candidate 
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FIG. 1: (a) Invariant mass distributions for the A'^ / can- 
didates in selected events, reconstructed in the pKtt (gray) 
and pK'^ (black) decay modes, (b) Invariant mass of the A~ 
candidate vs. that of the A'^ candidate in the same event, in 
5 MeV/c^ square bins. 



and at least one A~ candidate, requiring each candidate 
to have mass within 190 MeV/c^ of the fitted A'^ peak. 
We then form yl+ A^ pairs provided that they have no 
common tracks in their decay chains. For these 21,000 
pairs we show the candidate pK~iT'^ and pK^^ invari- 
ant mass distributions in Fig. [1^. Clear yl+ signals are 
visible over modest backgrounds. The peak mass val- 
ues, rates, and momentum distributions are consistent 
with previous measurements 12, 14, 1^. We plot the in- 
variant mass of the A~ candidate versus that of the A'^ 
candidate in Fig. [TJj. Horizontal and vertical bands are 
visible, corresponding to events with a real A~ or yl+, re- 
spectively, and there is a substantial enhancement where 
they overlap. 

The opening angle 9 between the A'^ and yl^ momenta 
in the cm. frame is sensitive to their production mech- 
anism. We expect A'^ A~ pairs from gluon splitting 
(e+e~— > qqg — > qqcc) or e+e^— ccg events with a very 
hard gluon to have relatively small 0, but also a sup- 
pressed selection efficiency due to the p* requirement. 
In the 21,000 events selected, 9 values are concentrated 



near 180°, consistent with dominance of 2-jet e+e"— >-cc 
events. Only seven events have 9 < 90°, one of which 
is in the signal region defined below. Since the small-6' 
background may have different characteristics from that 
at large 9, we require 6'>90°. 

About 3% of the events have two A'^ (or two A~) 
candidates, due to the two pK~ Tr~^ combinations in 
the decay chains 17++ -)■ yl+(pif ~7r+)7r+ and yl*+ ^• 
yl+(pif ~7r+)7r+7r^. We include all combinations in the 
sample and account for the kinematic overlap through the 
background subtraction. We define a circular A^A'^X 
signal region centered at our peak mass values with a 
radius of 12 MeV/c^, which contains 919 entries. Using 
the single- yl+/yl~ bands [13[, we estimate an expected 
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FIG. 2: Background-subtracted distributions for the 649 
A'^A~X events in the data: (a) the numbers of additional 
tracks, identified and identified p/p; and (b) missing mass, 
with imaginary masses given negative real values. Most events 
have no identified or p/p and the corresponding zero- 
multiplicity points are off the vertical scale in (a). 



estimate 25 ± 1 expected background events with fake yl+ 
and A~ , giving a A'^A~X signal of = 649 ± 35 

events. 

We can calculate an expected number of signal events. 



background in the signal region of 245 ± 5 events with 
one real yl+ or A~ and one fake. Using events with both 
masses at least 40 MeV/c^ from the fitted yl+ mass, we 



under the assumption that the c and c hadron types 
are uncorrelated so that all signal events are four-baryon 
events. Then Uexp = Cnl/ANcc, where ni = 420,000 
is the number of single A'^/A~ observed in the data, 
A^cc^SxlO* is the number of e+e^— >cc events expected 
for our integrated luminosity, and the factor C accounts 
for the correlation between the /!+ and A^ reconstruction 
efficiencies. This formulation is independent of the yl+ 
branching fractions and average efficiencies. In the simple 
case where the efficiencies of the yl+ and A~ in A'^A'^X 
events are uncorrelated, no correction is needed (C = 1) 
and Uexp = nf/4Ncc- More generally, < C < 1/e for an 
average acceptance times efficiency of e: in the extreme 
case of maximal correlation C — 1/e, and in the extreme 
case of maximal anticorrelation riexp = C = 0. At EaBAR 
there might be correlations because of the asymmetric 
beam energies and detector layout. We evaluate this cor- 
rection using the JETSET, HERWIG, and UCLA mod- 
els, adjusting their charm fragmentation parameters and 
reweighting the resulting p* distributions to reproduce 
our measured distribution for inclusive /!+ [l^ . Com- 
bined with smooth parametrizations of our efficiencies as 
functions of momentum and polar angle, the models give 
values of C ranging from 0.63 to 1.65, with a mean of 1.05. 
Even allowing for the large model dependence, the full 
range of Uexp — 100-250 events is well below the observed 
649 ± 35, confirming the enhanced rate N^+-^- /uexp ~ 4 
reported by the CLEO Collaboration [iq| . 

We investigate the structure of the A'^A'^X events 
using the 71+ and A~ candidates along with additional 
charged tracks that have at least 10 points measured in 



the DCH, 5 in the SVT, and extrapolate within 5 mm 
of the beam axis. We subtract appropriately scaled dis- 
tributions in the background regions from those in the 
signal region to obtain distributions for A'^A~X events. 
Figure [5^ shows the distribution of the number of ad- 
ditional tracks, as well as the numbers of identified 
and p/p among them. Were each c baryon compensated 
by a light antibaryon, then — assuming that half the an- 
tibaryons have an antiproton in the final state and ac- 
counting for p/p detection efficiency — we would expect 
45% of these events to contain one identified p/p and 
another 20% to contain both an identified p and a p; 
we observe only 3.4% and 0.6%, respectively. Figure [5]d 
shows the distribution of missing mass, calculated from 
the four- momenta of the initial e+ and e~, the recon- 
structed and A'^ , and all additional tracks interpreted 
as pions. A typical NcifiXnN c2 event, containing both a 
neutron and an antineutron, would have a missing mass 
well in excess of 2 GeV/c^ . 

The distributions in Fig. [5] indicate that the majority 
of the A'^A~X events do not contain additional baryons, 
and therefore that the conservation of baryon number 
is realized with the primary c and c hadrons. In the 
background-subtracted sample of 649 ± 35 AfA~X sig- 
nal events, there are 28 ± 6 additional identified p/p can- 
didates. These p/p candidates include background from 
two main sources: interactions in the detector material 
and misidentified pions or kaons. We expect 5 protons 
from material interactions. We also expect about 12 pi- 
ons or kaons misidentified as protons, based on the num- 
bers and momenta of the observed additional tt^ and 
tracks. In cross-checks these expectations are found 
to be consistent with the data within uncertainties: there 
are 8 ± 4 more identified p than p (with the excess at- 
tributed to material interations) , and there are 7 ± 3 
events seen with exactly one additional identified p/p 
and an event missing mass below 750 MeV/c^ (incon- 
sistent with a missing second baryon, and so attributed 
to a misidentified kaon or pion). Subtracting the ex- 
pected contributions from these two background sources, 
correcting for efficiency, and assuming equal p and n pro- 
duction rates, we estimate that we observe 13 ± 8 true 
four-baryon events. This is well below the rate of 100 
to 250 four-baryon events expected for uncorrelated pro- 
duction, let alone the observed rate of 649 ± 35 events, 
indicating that the four-baryon process is strongly sup- 
pressed and that the primary production process domi- 
nates. 

None of the reconstructed events is consistent with 
the two-body process e+e"— >■ yl+vl". However, the sig- 
nal could arise from the pair-production of c baryons if 
one or both are excited states that decay to A'^/A~: 
e+e~ ^ NciNc2^ A+A-X. Combining yl+/Z J candi- 
dates with one or two additional tracks assigned the pion 
mass hypothesis gives the invariant mass distributions 
in Fig. [3] The points represent sideband-subtracted sig- 
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FIG. 3: Invariant mass distributions for (a) A'^tv^ and A'^tv^ 
and (b) A'^ / A~ tt'^ -k~ combinations. The points with errors 
represent the background-subtracted A'^A^X events, and the 
weighted histograms are from the single-ylj/yl J sidebands. 



nal events and the histograms the single- yl+/yl J side- 
bands with entries reweighted to reproduce the num- 
ber of the yljf/yljT in signal events and their momentum 
and polar angle distributions in the lab frame. Peaks 
are visible in the sideband data for the 17^^^° (2455), 
i:++/°(2520), and the excited A+ states at 2593, 2625, 
2765 and 2880 MeV/c^. We find no unexpected peaks in 
our yl+7r(7r), A'^K or A^p mass distributions. The points 
are consistent with the histograms, indicating similar 
c baryon compositions in the two event types. Only two 
events are kinematically consistent with NdN c2 

for these known Nc- Distributions of 6 and the decay 
angles in the vI+tt rest frames are consistent with multi- 
hadron events, and not with very heavy states decaying 
into a yl+ and more than two pions. We conclude that 
e+e"— >■ NciNc2 processes represent a small fraction of 
our sample. From the fits in Fig. [3J we estimate that 
35 ± 3% of the A+ and 29 ± 2% of the additional pions 
in our sample are decay products of heavier c baryons. 

Having established the presence of a category of events 
containing a c baryon, a c baryon, no other (anti)baryons, 
and several intermediate mesons, we study the num- 
ber and structure of these mesons. We exclude events 
with an identified p/p or a missing-mass-squared below 
—0.25 GeV^/c^. We estimate that the sample contains 
a further 5 ± 5 four-baryon events in which no p/p is 
detected; we take these to have the same distributions 
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as the events with an identified p/p and subtract an ap- 
propriately scaled contribution to correct for them. In 
this sample of 619 ± 35 events we study a number of 
quantities including the / A~ and additional track mo- 
menta, polar angles, rapidities and opening angles. Their 
inclusive distributions are quite similar to those in the 
single- yl+/yl^ sample and similar to those in all hadronic 
events. In particular, signing the thrust axis such that 
the rapidity is positive, the and A'^ rapidities 
cluster near -1-1.1 and —1.1 units, respectively, with the 
additional tracks of each charge distributed broadly and 
symmetrically in between. 

These 619 events contain only 45 ± 10 identified 
of which about 20 are expected to be misidentified pious. 
The events show no mass peak for candidates recon- 
structed from pairs of tracks not included in the or 
A~ (including tracks that do not extrapolate within 5 mm 
of the beam axis). The K'.i: ratio is thus much lower than 
the value 0.3 typical of hadronic events, which might be 
due to the limited energy available and the fact that our 
c baryons are non-strange (the lighter c-s baryons do not 
decay into A'^). The tt+tt", fC^Tr^, and K^K~ invariant 
mass distributions show no significant resonant structure; 
in particular there is no evidence for the . This implies 
a vector:pseudoscalar meson ratio much lower than the 
value near 1 typical of hadronic events, and suggests that 
most tracks not from c baryon decays represent distinct 
intermediate mesons. 

The intermediate meson multiplicity is distributed 
broadly. We verify that the contribution from decays of 
heavier c baryons is not concentrated in any particular 
region in Fig. [2^, but due to the limited sample size we do 
not attempt to correct the distribution. We observe an 
average of 2.7 additional charged tracks per event. Cor- 
recting for c baryon decays and tracking efficiency gives 
2.6 ± 0.2 charged intermediate mesons per event, where 
the uncertainty includes both statistical and systematic 
effects. The uncertainty is dominated by the track accep- 
tance in these events, evaluated with a set of simulations 
based on the observed tt* and distributions. On av- 
erage, the c and c baryons carry 75% of the event energy, 
and the intermediate charged mesons account for about 
65% of the remainder. This and the broad distribution 
of missing masses in Fig. [2)d suggest the presence of ad- 
ditional neutral mesons. If intermediate 7r° are produced 
at half the tt^ rate, as in typical hadronic events, the av- 
erage intermediate meson multiplicity would be 3.9 ±0.3. 

The new type of event observed in our data might be 
explained by either primary diquark-antidiquark produc- 
tion or the production of multiple intermediate mesons 
between a baryon and antibaryon. Neither the JETSET 
nor the HERWIG model produces events of the type ob- 
served, although both might be adapted to include one 
or both of the above processes. JETSET does produce 
NciMN c2 events, where M is a single meson, often a vec- 
tor decaying into two or three pions, but the event char- 



acteristics are far from consistent with the data. Multi- 
ple intermediate meson processes occur naturally in the 
UCLA model, which also predicts an enhanced A'^A~X 
fraction due to events of this type, with suppressions 
of kaons and vector mesons. The version of the UCLA 
model used does not describe the observed events in de- 
tail, having an average of only 1.8 intermediate mesons 
with a distribution peaked at low values, but the results 
presented here should encourage development of this and 
other relevant models. 

In summary, we isolate a sample of 649 ±35 e"'"e~— >cc 
events containing both a A'^ and a A'^ with high momen- 
tum in opposite hemispheres, and study these events in 
detail. The number of events is estimated to be about 4 
times that expected if the leading c and c hadron types 
are uncorrelated, confirming an observation by the CLEO 
Collaboration. Taking advantage of the particle identi- 
fication capabilities of the BaBAR detector and the large 
data sample, we are further able to establish that al- 
most all of these events contain no additional baryons. 
They do contain 2.6 ± 0.2 additional charged intermedi- 
ate mesons on average, and events with zero additional 
mesons do not contribute significantly. Our event sample 
exhibits distributions of momentum, angle, rapidity, and 
c baryon type similar to those in typical hadronic events, 
but contains fewer kaons and vector mesons. This is di- 
rect evidence for a new class of multihadron events, in 
which baryon number is conserved by a leading baryon 
and antibaryon, rather than locally along the hadroniza- 
tion chain. 
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